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Abstract 

The decay B —* /o(980)-K" is studied within the framework of QCD factorization and the two- 
quark scenario for /o(980). There are two distinct penguin contributions and their interference 
depends on the unknown mixing angle of strange and nonstrange quark contents of /o(980): 
destructive for < 6 < tt/2 and constructive for tt/2 < 9 < tt. The QCD sum rule method 
is applied to evaluate the leading-twist light-cone distribution amplitudes and the scalar decay 
constant of /o- We conclude that the short-distance approach is not adequate to explain the 
observed large rates of foK~ and fo~K°. Among many possible subleading corrections, we study 
and estimate the contributions from the three-parton Fock states of the /o and from the intrinsic 
gluon inside the B meson. It is found that the spectator gluon of the B meson may play an eminent 
role for the enhancement of /o(980)_RT . We point out that if /o(980) is a four-quark state as widely 
perceived, there will exist extra diagrams contributing to B — » Jq(980)K. However, in practice it is 
difficult to make quantitative predictions based on the four-quark picture for /o(980) as it involves 
additional nonfactorizable contributions that are difficult to estimate and the calculations of the 
decay constant and form factors of /o(980) are beyond the conventional quark model. 



1 



I. INTRODUCTION 



The decay of the B meson into a scalar meson /o(980) was first measured by Belle [l| in the 
charged B decays to K^it^tt^ and a large branching fraction product for the fo(980)K ± final states 
was found. A recent updated result by Belle yields 

B(B + -» / (980)K+ -> ir + ir-K + ) = (7.55 ± 1.24±£;Jjj) x 10~ 6 . (1.1) 

The Belle result is subsequently confirmed by the BaBar measurement [3J: 

B{B + -> /o(980)ET + -> ^ + vr-K+) = (9.2 ± 1.2±|;^) x 10~ 6 . (1.2) 

A recent BaBar analysis of B^ — > i^^vr^vr 111 gives a very similar result: (9.2 ± 1.5 ± 0.8) x 10~ 6 0|. 
The world average is then given by 

B(B + -» / (980)K + -»■ ^+7^^+ ) = (8.49±1;||) x 10~ 6 . (1.3) 

BaBar has also measured the neutral mode B° — ► /o(980)iT° with the result 0| 

i3(5° -> / (980)K° -> ^ + ^if°) = (6.0 ± 0.9 ± 1.3) x 10" 6 . (1.4) 

This channel is of special interest as possible indications of New Physics beyond the Standard 
Model (SM) may be observed in the time-dependent CP asymmetries in the penguin-dominated 
B decays such as B° -> f (980)K s . The mixing- induced CP-violation parameter S is expected to 
be very close to — sin (5 in the SM. The most recent measurements by BaBar and Belle yield 



sin/ 3 (/o^) = (°- 95 -- ±0 - 10 BaBar J 
HyJ ' 1-0.47 ±0.41 ±0.08 Belle |g. 



The deviation from sin 2/3 = 0.726 ± 0.037 |9j measured from B — > J/ipKs may hint at a possible 
New Physics. 

The absolute branching ratios for B — ► JqK depends critically on the branching fraction of 
/o(980) — ► 7T7T. For this purpose, we use the results from the most recent analysis of lOj, namely, 
IV* = 64 ± 8 MeV and r tot = 80 ± 10 MeV for / (980), to obtain B(/ (980) -» ttvt) = 0.80 ± 0.14. 
1 This leads to 

B(B+ -> /o(980)i^ + ) » (15.9±g;|) x HT 6 , 

i3(S° -» /o(980)i^°) « (11.3 ± 3.6) x HT 6 . (1.6) 

Comparing with the averag ed branching ratios, (12.1 ± 0.8) x 10~ 6 for B+ tt°K + and (11.5 ± 
1.0) x 10~ 6 for B° -> vr°K fl, we see that / (980)is: + > tt°K + and M980)i^° ra tt K . 



Theoretically, the decay B — ► /o-fT has been studied in 



121 and 131 within the framework of 



the pQCD approach based on the hr factorization theorem. It is found that the branching ratio is 



1 The ratio R = B(fo -> Tr+ir-)/B(f -> K+K~) w 7.1 is consistent with the result of i? > 2.6± : 
inferred from the Belle measurements of B(B + — > /o(980)iir + — > 7r + 7r~A" + ) [see Eq. Ijl.ip ] and B(B + 
fo{98Q)K+ -> K+K-K+) < 2.9 x 10" 6 Q. 
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of order 5 x 10~ 6 (see Fig. 2 of ["U^b which is smaller than the measured value of B(B + — ► /qK + ) 
by a factor of 3. In the present paper we shall re-examine this decay within the QCD factorization 
approach 1J, Il5l . Il€ ] . 



The underlying structure of the parity-even meson /o(980) is still controversial and not clear: It 
can be a conventional 2-quark P-wave state or a S-wave meson made of four quarks qqqq. It will 
be interesting to see if these two different scenarios for /o(980) can be tested in B — > /qK decays, 
an issue which we will address briefly in this work. 

There are two distinct penguin contributions to B — > foK: one is related to the u quark 
component of /o an d the other to the strange quark content of /o- Owing to the cancellation 
between two penguin terms, the former penguin contribution is severely suppressed, while the 
latter depends on the unknown scalar decay constant of Jq. Based on the sum rule approach, 
we shall show that the magnitude of the fo scalar decay constant is sunstantially larger than the 
corresponding decay constant of pseudoscalar mesons and hence the f$K rate can be comparable 
to the tt°K one. However, the predicted branching ratio of B — > f$K is still lower than experiment 
by ~ 45% for the f$K~ mode and ~ 30% for foK . There are several possible mechanisms such 
as final-state interactions, flavor-singlet contributions, large annihilation contributions that may 
account for the enhancement of f$K. In the present work, we will focus on the subleading effects 
arising from the three-parton Fock states of the /o and from the spectator gluon inside the B meson. 

It is commonly assumed that only the valence quarks of the initial and final state hadrons 
participate in the decays. Nevertheless, a real hadron in QCD language should be described by a 
set of Fock states for which each state has the same quantum number as the hadron. For instance, 

IB') = *pg\bu) + ^ g \bug) + ipgjbuqq) + ^Jbucc) + . . . , 

l/o) = Vnn\nn) + tpig\ss) + ^l° ng \nng) + ipi§ g \ssg) + XX™*»?d nf ^ + ip%q\ssqq)) + ... , 

= ^fu\su) + ip^ g \sug) + ^fuqq\suqq) + ... , (1.7) 

where nn = (uu + dd)/\/2. The extra gluon(s) or quark pair(s) appearing in higher Fock states 
are the results of QCD interactions. The cc pairs due to a single gluon splitting g — > cc, described 
by the Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution equation, 2 are basically extrinsic to the 
bound-state nature of the hadron. In contrast, the cc pairs, which are entangled through multiply 
gluonic interactions to the valence quarks, at least via g*g* — > cc, are intrinsic to the hadronic 
structure EJ. It has been estimated that the intrinsic charm probability in the proton is 

?S 1% |3 The intrinsic charm may explain the pir puzzle in J/tp(ip') decays (3]. 

The study of the intrinsic charm component in the B meson is interesting and has attracted a 
great deal of attentions 1£ , 2(3, 21 ] . Although the extrinsic charm pairs carry only small momentum 



2 We should remind the readers that the light-cone wave functions in Eq. (|1 .T|) are Lorentz invariant using 
the light-cone quantization in the light-cone gauge A + = 17]. However, here we adopt the equal-time 
quantization and choose the B rest frame. Each Fock-state wave function of the B meson is no longer boost 
invariant. In contrast, the energetic /o and K produced in B decays are represented in terms of light-cone 
distribution amplitudes (LCDAs) in conformal expansion Q|. Strictly speaking, the GLAP equation can 
be applied only to the partons in the infinite momentum frame or to the light-cone quantization framework. 
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fraction of hadrons, it has been argued that the intrinsic cc pair could share 20% momentum fraction 
of the parent B meson (20I ] . As discussed in , the intrinsic charm may give rise to 20% corrections 
to B — > Kit amplitudes. Consistently, \^buccl ^bu? integrating over the phase space is estimated to 
be roughly of order (20%) 2 = 4% [21| . compared to the charm content < 1% in the proton. Similar 
to the intrinsic charm case, one may ask how large is the intrinsic gluon within the B meson. The 
gluon content of the B meson is intimately related to the energy of light degrees of freedom (or the 
so-called "brown muck"), namely, A = tub — rrib- Even though A > 450 MeV is larger than the 
constituent light quark mass, it may suggest that the intrinsic gluon within the B meson cannot 
be neglected. To estimate the possible contributions originated from the intrinsic gluon within the 
B meson, we introduce the three-particle operator by a q„G nf i,u, which can couple to the B meson. 



Similar definitions for it and p can be found in the literature |2J, |25j. Hence, the intrinsic gluon 
effects can be estimated. The detailed calculation is shown in Sec. Ill CI As a result, we find that 
the intrinsic gluon may give rise to 20-40% corrections to the decay amplitudes. In contrast, we 
shall show in Sec. Ill Bl that the extrinsic gluon effect due to the splitting b —* bg is negligible, which 
is described by the GLAP evolution equation in the infinite momentum frame. 

The layout of the present paper is organized as follows. In Sec. II we shall first study B — * foK 
decays within the framework of QCD factorization and discuss the internal structure of the /o(980). 
We then continue to explore the subleading effects arising from the three-parton Fock states of the 
/o and from the spectator gluon inside the B meson. We present numerical results and discussions 
in Sec. Ill and give conclusions in Sec. IV. Appendices A and B are devoted to the determination 
of the scalar decay constant and distribution amplitudes of /o, respectively, Appendix C outlines 

the sum-rule calculation for the decay constant 5 2 fq and Appendix D for the intrinsic gluon content 

Jo 

of the B meson. 

II. B -> /o (980) AT DECAY AMPLITUDES 

A. B — > f (980)K decay amplitudes in QCD factorization of leading Fock states 

1. Framework 

To proceed we first discuss the decay constants and form factors. The decay constants are 
defined by 

(K(p)\A^\0) = -if KPlM </ |Vg0>=0, (fo\qq\0)=m f J q . (2.1) 

The scalar meson /o(980) cannot be produced via the vector current owing to charge conjugation 
invariance or conservation of vector current. The scalar decay constant f q will be discussed later. 
Form factors for B — > P and B — > S transitions (P: pseudoscalar meson, S: scalar meson) are 



defined by |2C ] 



(P(PP)W,\B( PB )) = (pB,+PP,-^^ q \ F BPtf ) + ?^& qiiF BPtf 



(2.2) 
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(a) 




(b) 



FIG. 1: Penguin contributions to B — ► fo(980)K in the 2-quark picture for /q(980). 



where = (pb - Pp)/i, and |23] 3 
(5(ps)|^|B(pb)> 



+- 



m B~ m S„ TrBSrJl 



(2.3) 



The penguin-dominated B~ — * foK~ receive two differen t ty pes of penguin contributions as 
depicted in Fig. 1. Within the framework of QCD factorization |l4j |. the B — > JqK decay amplitudes 
read 4 

A(B-^f K-) = ~ E A 4 [ai(/o^-)^ + K-r^)(/ K-) + K -r^)(/o^-) 



A(B° -> /o^°) 



(m^)(m| - mf ) + (2ag - ag)(Jir/o)/. ^fif*(m$,)(mf, - ml 

E /B/i</ 9 k(/0^-) + ^ EW (/0^" 
q=u,s 

-% E V 

V ^ p=u,c 



K J 



(aS - ffag)(/oX°) - J« " ifaS)(/aK°) 



E fsflifq 



q=u,s 



(2.4) 



where X p = V p bV* s , r^(fi) = 2m 2 K /[rrib(lj){rn u (ii) + m s (/j,))], and weak annihilation contributions 
described by 63 and 63EW terms will be discussed shortly. In Eq. (|2.4|) the superscript u of the 
form factor F Q reminds us that it is the u quark component of /q involved in the form factor 
transition [see Fig. 1(a)]. In contrast, the subscript s of the decay constant f s indicates that it is 
the strange quark content of fo responsible for the penguin contribution of Fig. 1(b). 



3 As shown in [2j|, a factor of (— i) is needed in Eq. (|2.3|l in order for the B — ► S form factors to be positive. 
This also can be checked from heavy quark symmetry [27|. 

4 The relative sign between ag(foK)fK and a^(Kfo)f s terms in the decay amplitude of B — > /qK obtained 
in [3| is opposite to ours. 
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The effective parameters a? with p = u, c in Eq. (|2.4|) can be calculated in the QCD factoriza- 
tion approach 14 1. They are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions Q, Q] 



af(M x M 2 ) 



. Ci±l , Q±l CFOi s f 



AL 



TV,. 4vr 



4tt 



Vi(M 2 ) + —H i {M 1 M 2 ) + Pf (M 2 ) 



2V, 



(2.5) 



where i = 10, the upper (lower) signs apply when i is odd (even), c, are the Wilson coefficients, 

C F = (N 2 - l)/(2N e ) with iV c = 3, M 2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi(M 2 ) account for vertex corrections, i!j(MiM 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and PjjM?.) for penguin contractions. The explicit expressions of these 



quantities can be found in I4I 

fufB 



1( 



The hard spectator function H reads 



Jo P jo 4 Jo V rn b £ J0 



1 dfy 



F S/o "(0)m 2 



for i = 1,4, 10 and fT, 



for i 



(2.6) 



6, 8. where £ = 1 — £ and rj = 1 — n. As for the parameters 



Og g(Ar°/o) appearing in Eq. ([2.4)1 . they have the same expressions as agg(/o^°) except that the 
penguin function Gk (see Eq. (55) of 14j) is replaced by Gf and by $j q . 

Weak annihilation contributions to B — > foK are described by the terms 63 and 63EW m Eq. 
Q2.4|) which have the expressions 



b q - ^ 
63 " AT| 



C3 ^r+c 5 (4 



°3EW 



Cp 



C9 4 q) + c 7 (4 q) + 4 iq) : 



) + N c c e A f 3 (q) 



(2.7) 



where A 3 is the factorizable annihilation amplitude induced from {S — P){S + P) operator and A\ 3 
are nonfactorizable ones induced from (V — A)(V — A) and (S— P)(S+P) operators, respectively. It 
is evident that the dominant annihilation contribution arises from the factorizable penguin-induced 



I— t 



annihilation characterized by A 3 . Their explicit expressions are given by (see also |l6j ] 

1 1 



ttq s / dxdy\$ K (x)$ { ^{y) 



y(l - xy) 



+ 



x 2 y 



2r*-^<S> p K {x)$ 



A 



A 



f{i) 



7ra, 



ira. 



dxdy 



dxdy 



2m 



fa 



{ m b 
f 2m fo 



$k(x)$> 



fo ( 



2j/ 



x mb 



fo 



2 

xy 



<^ K (x)^ p (y) 



xy(l - xy) 
2(1 + x) _ 
x 2 y 



+ r^>{y)$ p K {x 



2x 



I nib 11 v / Jo w x 2y x ju - .... X y± 

where q = u, s, x = 1 — x, y = 1 — y, <&m (&m) is the twist- 2 (twist-3) light-cone distribution 
amplitude of the meson M. 

Although the parameters cij(z 7^ 6, 8) and a^^r x are formally renormalization scale and 75 scheme 
independent, in practice there exists some residual scale dependence in ctj(/i) to finite order. To be 
specific, we shall evaluate the vertex corrections to the decay amplitude at the scale \i = nib/2. In 
contrast, as stressed in Q|, the hard spectator and annihilation contributions should be evaluated 



xy(l - xy) 
2(1 +y) 



(2. 
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at the hard-collinear scale ph = vT^h with ~ 500 MeV. There is one more serious complication 
about these contributions; that is, while QCD factorization predictions are model independent 
in the rrif, — > oo limit, power corrections always involve troublesome endpoint divergences. For 
example, the annihilation amplitude has endpoint divergences even at twist-2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and posses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow 14] to parametrize the endpoint divergence Xa = Jq dx/x in the annihilation 
diagram as 

X A =\n(^yi + p A e i ^), (2.9) 

where pa is a complex parameter < < 1. Likewise, the endpoint divergence Xh in the hard 
spectator contributions can be parametrized in a similar way. 

Note that 04 and a% penguin terms contribute constructively to tt°K~ but destructively to 
foK~~. Therefore, the contribution to B — > foK from Fig. 1(a) will be severely suppressed. The 
contribution from Fig. 1(b) is suppressed by m/ /m&. Hence, it is naively expected that the foK 
rate is much smaller than the ir®K one. However, as we shall see below, the scale dependent decay 
constant f s is much larger than f n . As a consequence, the branching ratio of B — > foK turns out 
to be comparable to B — > 7r°K. 



2. Quark structure of fo(980) and input parameters 

It is known that the underlying structure of scalar mesons is not well established theoretically 
(for a review, see e.g. 29, 3(3, 31 ]). It has been suggested that the light scalars below or near 1 GeV- 
the isoscalars /o(600) (or a), /o(980), the isodoublet 1^(800) (or k) and the isovector ao(980)-form 
an SU(3) flavor nonet, while scalar mesons above 1 GeV, na mely , /o(1370), ao(1450), ^(1430) and 
/o(1500)//o(1710), form another nonet. A consistent picture |3J| provided by the data suggests that 
the scalar meson states above 1 GeV can be identified as a conventional qq nonet with some possible 
glu e content, whereas the light scalar mesons below or near 1 GeV form predominately a qqqq nonet 



32 



331 ] with a possible mixing with + qq and glueball states. This is understandable because in 
the qq quark model, the + meson has a unit of orbital angular momentum and hence it should 
have a higher mass above 1 GeV. On the contrary, four quarks q 2 q 2 can form a + meson without 
introducing a unit of orbital angular momentum. Moreover, color and spin dependent interactions 
favor a flavor nonet configuration with attraction between the qq and qq pairs. Therefore, the + 
q 2 q 2 nonet has a mass near or below 1 GeV. This four-quark scenario explains naturally the mass 
degeneracy of /o(980) and ao(980), the broader decay widths of <r(600) and k(800) than /o(980) 
and ao(980), and the large coupling of /o(980) and ao(980) to KK. 

While the above-mentioned four-quark assignment of /o(980) is certainly plausible when the light 
scalar meson is produced in low-energy reactions, one may wonder if the energetic /o(980) produced 
in B decays is dominated by the four-quark configuration as it requires to pick up two energetic 
quark-antiquark pairs to form a fast-moving light four-quark scalar meson. The Fock states of 
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B 



B 




u 



(b) 

FIG. 2: Penguin contributions to B~ — > /o(980)i^~in the 4-quark picture for /o(980). 

/o(980) consists of qq, q 2 q 2 , qqg etc. [see Eq. (|1.7|) ]. It is thus expected that the distribution 
amplitude of /o would be smaller in the four-quark model than in the two-quark picture. Then one 
will not be able to explain the observed B — > fo(980)K decays. 

Nevertheless, as pointed out in [3^], the number of the quark diagrams for the penguin contri- 
butions to B — > fo(98Q)K (see Fig. |2J) in the four-quark scheme for /o(980) is two times as many 
as that in the usual 2-quark picture (see Fig. That is, besides the factorizable diagrams in 
Fig. Eta), there exist two more nonfactorizable contributions depicted in Fig. |2^b). Therefore, a 
priori there is no reason that the B — > /o(980)-K" rate will be suppressed if fo is a four-quark state. 
However, in practice, it is difficult to give quantitative predictions based on this scenario as the 
nonfactorizable diagrams are usually not amenable. Moreover, even for the factorizable contribu- 
tions, the calculation of the /o(980) decay constant and its form factors is beyond the conventional 
quark model, though attempt has been made in In order to make quantitative calculations for 
B — > fo(980)K, we shall assume the conventional 2-quark description of the light scalar mesons. 

In the naive 2-quark picture, /o(980) is purely an ss state and this is supported by the data of 
Df — > /ovr + and (ft — ► /o7 implying the copious /o(980) production via its ss component. However, 
there also exist some experimental evidences indicating that /o(980) is not purely an ss state. First, 
the observation of T(J '/if; — > Jquj) k, ^T(J/i(j — > focft) |35( clearly indicates the existence of the non- 
strange and strange quark content in /o(980). Second, the fact that /o(980) and ao(980) have 
similar widths and that the /o width is dominated by irir also suggests the composition of uu and 
dd pairs in /o(980); that is, /o(980) — > irn should not be OZI suppressed relative to ao(980) — ► nr]. 
Therefore, isoscalars <t(600) and /o must have a mixing 

|/o (980)) = \ss) cos 9 + I nn) sin 0, M600)) = -\ss) sinfl + \nn) cos 9, (2.10) 

with nn = (uu + dd) / \[2. The distribution amplitudes $ s = $^ and <I> n = 3>^P' corresponding to 
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/o = ss and f$ = nn, respectively, are 

</ »|^h^(0)|0> = pj n ^ dxe^ z ^ n (x), 

Jo 

(f^pMz) lfl s(0)\0) = P J S f 1 dxe ix ^ s (x), 

Jo 

</?(p)|n(z)n(0)|0) = m^f n £dxe ix ^ n (x), 

(m P Mz)s(0)\0) = mff s j\xe^ z ^ s {x\ (2.11) 

with f q being denned by 

(f^\qq\0)=m^f q . (2.12) 

They satisfy the relations <& n ,s( x ) = ~~ ^n,s(l ~ x ) due to charge conjugation invariance (that is, the 
distribution amplitude vanishes at x = 1/2) and <J?^ )S (x) = $^, s (l — x) so that dx Q ntS (x) = 0. 
For the scalar meson made of qq, its general distribution amplitude has the form 



$s(x,/i) = &c(l 



B (v)+ B m {n)CU 2 (2x-l) 



m=l 



(2.13) 



3 /2 

where -B m are Gegenbauer moments and Cm are the Gegenbauer polynomials. For the isosinglet 
scalar mesons a and /o, Bq = and only the odd Gegenbauer polynomials contribute. Hence, the 
light-cone distribution amplitudes (LCDAs) for /o read 

<S> n>s (x,tM) = 6x(l-x) ]T Bt s) (^)C 3 J\2x-l). (2.14) 

m=l,3,5,--- 

Using the QCD sum rule technique we have evaluated in Appendix B the Gegenbauer moments 
up to m = 5. In particular, we obtain B[ n) = -0.92 ± 0.08 for $ n and b[ s) = 0.8B[ n) for $ s at 
^ = 1 GeV. Our results are consistent with the empirical result of \Bi\ = 1.1 inferred from the 
analysis in [3^. As we shall see below, the negative sign of B\ is indeed strongly favored when 
the predicted foK rates are confronted with the data. As for the twist-3 distribution amplitude 
$j o (x), its asymptotic form is the same as the light pseudoscalar meson to the leading conformal 
expansion |l8(. Hence, we take 

K,.( x ) = L ( 2 - 15 ) 

The asymptotic forms for kaon twist-2 and twist-3 distribution amplitudes are 

$k(x) = 6x(l - x), $ p K (x) = 1. (2.16) 
In the qq description of /o(980), it follows from that 

D— -£ ! R — fUU r-,n r I r0 fdd 

F b fo = -j=^6 Fq f ° , Fq f0 = ^sin9 F B f ° , (2.17) 

where the superscript qq denotes the quark content of fo involved in the transition. The form 
'actor for B to the scalar meson transition has been calculated in the covariant light-front model 
3. From Table VI of 



Bf qq 

22J, it is clear that F (0) with qq = uu or dd is of order 0.25 which is 
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very similar to F^*(0). As we shall see below, a precise estimate of the B — > /o(980) form factor 
is not important because the contribution from Fig. ^a) is largely suppressed owing to the large 
cancellation between the 04 and a& penguin terms. Based on the QCD sum-rule method, the scalar 



decay constant f s defined in Eq. (|2.12|) has been estimated in |37( and |38(] with similar results, 
namely, f s ~ 0.18 GeV at a typical hadronic scale. Notice that, in contrast to the pseudoscalar 
meson decay constants, the scalar decay constant f q is scale dependent. Taking into account the 
scale dependence of f q and radiative corrections to the quark loops in the OPE series, in Appendix 
A we have made a careful evaluation of the scalar decay constant using the sum rule approach 
and found / s (/i = 1 GeV) = 0.33 GeV and / s (/i = 2.1 GeV) = 0.39 GeV [see Eqs. (ESJ) andln the 
two-quark scenario for /o(980), the decay constants f n>s are related to f n>s defined in Eq. (|2.1|) via 
(|A10|) ] and similar results for f n . 

(s) (n) 

Is = — -fs cos9, f n = -J^f nS m6, (2.18) 

m fo m fo 

where use has been made of Eqs. (|2.1|) . (|2.1U|) and H2.12j) . 

Experimental implications for the fo—cr mixing angle have been discussed in detail in [3^: 

J/ip -> f <f>, f u [39] = (34 ± 6)° or 9 = (146 ± 6)°, 

R = 4.03 ±0.14 [39] => 9 = (25.1 ± 0.5)° or 9 = (164.3 ± 0.2)°, 
R = 1.63 ± 0.46 [39] 9 = (42.3±|;f )° or 9 = (158 ± 2)°, 

<t> -> /07. fo 77 [40] => = (5±5)° or 6 = (138±6)°, 
QCD sum rules and / data [41] 6» = (27 ± 13)° or = (153 ± 13)°, 

QCD sum rules and a data [41] =^ 9 = (41 ± 11)° or = (139 ± 11)°, (2.19) 

where i? = g"j QK + K - / / S'j 07r + 7r - measures the ratio of the /o(980) coupling to K + K~ and 7r + 7r~. In 
short, 6 lies in the ranges of 25° < 9 < 40° and 140° < 6 < 165°. Note that the phenomenological 
analysis of the radiative decays 4> — > /o (980)7 an d /o(980) — > 77 favors the second solution, namely, 
9 = (138±6)°. 5 

B. Subleading QCD factorization amplitudes arising from 3-parton Fock states of 
the final state mesons 

We shall see in Sec. Ill that the leading QCD factorization contributions to B — * foK are not 
adequate to explain the observed large rate of foK~ and foK° rates. This leads us to contemplate 
some possible mechanisms for enhancement. In this subsection we study the subleading QCD 
factorization amplitudes arising from three-parton Fock states of the fo or kaon. As will be shown 
later, these corrections are of order 1/m^ and turn out to be negligible. The relevant diagrams are 
depicted in Fig. |21 The study of these corrections is motivated by the following observations: (i) As 







ight s 
' 42]. 



In the four-quark scenario for li ght scalar mesons, one can also define a similar f — a mixing angle. It 
has been shown that = 174.6° 
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B — 




B * 



(b-1) 




(a-2) 




(b-2) 



FIG. 3: The subleading QCD factorization amplitudes originating from the qqg Fock states of the 
final state mesons in the B — > f$K decays. 

shown in 43] , Fig. |2tb) could give significant corrections to B —* ujK and urn decays because of the 
enhancement by the softer gluon effect accompanied by the relatively large Wilson coefficients C4 
and cq. (ii) Consider Fig.|2ta) where the emitted meson is a pseudoscalar 14| or a vector meson |43j. 
Its total amplitude vanishes due to the mismatch of the G-parity between the emitted meson and 
the relevant 3-parton operators. In contrast, the G-parities of 3-parton operators do match with 
the fo meson. Moreover, the contributions can be further enhanced by the large Wilson coefficients 
C4 and cq. 

In the following calculation, we adopt the conventions D a = <9 a +ig s T"M^, G a p = (\j2)e a ^ lLV G llv 
and e 0123 = — 1. We also use the Fock-Schwinger gauge to ex pre ss the gluon field A a in terms of 



G" for ensuring the gauge-invariant nature of the results [43 , 44] 



-l 



dv vG® v (vx)x l 



The three-parton LCDAs of the fo meson are defined by 

(fo{p)\q( x )g s G fJ , u (vx)o- aj3 q(0)\0) 

Vvg m ) - p a (Pfiguf3 - PuQup)] 



-fl h \pp{p^gv 



(2.20) 



(fo(p)\q(x)^g s G a p(vx)q(0)\0) 

-ifq *>■''• - 



pp [ g^ 



Pa gpn 



xpp^ 



px 



px / \ px 

(p a xp- Pl 3X a ) f Vacf>l( ai )e^ aM ^ 



{foip)\q{x)^^g s G al3 (vx)q{Q)\&) = f q (p a 9/3n ~ Pp9a». 
x I Pa^e^( a »+" a ») - f q ^(p a xp-ppx ) [ Va'^ q 



+ 



1<1 \ e ipx{a u +va g ) 



(2.21) 



(2.22) 
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where T>a = da q da q da g 5{\ — a q — a q — a g ) with a q , a q , a g being the fractions of the fo momentum 
carried by the (/-quark, g-quark and gluon, respectively. Here 0^ is a twist-3 LCDA, and 0j_, </>jj, 
</>l 0| are twist-4 ones given by 



<f>lf (<*i) 



2 

360a q a q a g 



1 + w i,o 2 ( 7a g - 3 ) + w 2,o(2 - 4a q a q - 8a g + 8a g ) 



+ w ljl (3a; ? ag — 2a g + 3a 2 ,) 



0f(at) 



305 2 >ja 2 (l — a 9 ) 



-1205 2 q a q a q a g 



- + 2e(l - 2a 5 ) 



1 



305^9 (q 9 — a q )a 2 



+ e(l — 3a g 



+ 2e(l - 2a 9 



120(5j 9 e(a g — a q )a q a q a g 



in the conformal expansion (for a further investigation, see [1 
the case of the pion, and 



(2.23) 

We obtain e ~ 0.30, similar to 



5 2 fs (n = 2.1 GeV) = (0.08 ± 0.01) GeV 2 , 
6}g{fi = 2.1 GeV) = (0.09 ± 0.01) GeV 2 , 



(2.24) 



for qq = ss or qq = nn. A detailed calculation of 5 2 fq is shown in Appendix Consider the 
four-quark operator s a 7^(1 — 75)6/3 ?/37^(l ^75)^0 i n Fig- El where a and /3 are color indices. The 
result of Fig. I^a.l) is found to be 

d 4 k 



Fig.Ha.1) 



1 



dv v I d x 



{k\m^9sil =F J 5 )G au (vx)q(0)\0)x» 



x (^18(0)7*^(1 " l 5 )e l ^- k H(0)\B) 

= 2 Mk\s j^(1 - 75)6|5) / Do^-g-[M + 2#_) ± 2(Sj - 201)], 



(2.25) 



where all the components of the coordinate x should be taken into account in the calculation before 
the collinear approximation is applied. Hence, in Eq. (|2.25|1 . the exponent arising from G au (vx) 
is actually e ik3 ' xv , where k 9 is the gluon's momentum, and the resultant calculation can be easily 
performed in the momentum space with the substitution of x u — ► —(i/v)(d/dk 9 l ). Following the 
same lines, the result for Fig. |3Ia.2) reads 

d 4 k 



Fig.lHta.2) = -1 j da a J d 4 x J 

x (ir-|s(0)7 M (l-75) 



(2vr 4 ) 



(/ok"(0) 7 ^ s (l T 75)G Q ,(aa;)( ? (0)|0)x 1 ' 



A; 2 



(2.26) 



H 
3 



9 j* /o (l - %)b\B~) j Da^[W + 201) =F 2(^ - 2$.)], 
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where we have used the following trick 



fa* 2 „ 



2i7_ 







(2.27) 



p h x dp B - 

during the course of calculation and introduced two light-like vectors: (nL = 0), parallel to 
the momentum of fo, and (n+ = 0,n + • n_ = 2), so that p^ Q ~ (j)f ■ n + )n^_ = pf 0+ n^i, x 2 ~ 
x-x+,p f0 -7 =p /o+ 7_. 

In Fig.Efb), the emitted gluon becomes a parton of the kaon. To proceed, we first take G^ u (vx) ~ 
G fJiU (0)e wk 3 ' x and then adopt the collinear approximation k^ = {a g )pK in the final stage, where 
(a g ) is the averaged fraction of the kaon's momentum carried by the gluon. The calculation is 
straightforward and leads to 



Fig. EN 



1 



4AL 







dv 
Tr 



du(V2f n $ n {u) + f s $ s (u)){K-\d^(l - l5 )g s G uP b\B- 
M u (v 



(vk£ + up h ) 2 



(vkf + up h ) 2 



0, 



(2.28) 



where we have used the fact that & n (u) = — <3? n (ii), $> s (u) = —& s (u) with u = 1 — u. 

In short, the perturbative contributions coming from higher Fock states of the final state mesons 
to the decay amplitudes are (in units of Gf/V%) 



A(B~ - f K-) Figm 



A(B° -> f K% iM 



~(m% - m 2 K )F BK (m%) f ^-^-(^6^-4(2^1 + ^) 



V2 V 



(c 4 + c 6 



)(V2/„(20i + 0f) + / s (2^ 



+ 



+i(c 8 + c 10 )(^(2^ + - J s (20i + 0f) 



(2.29) 



However, the above result is numerically negligible. Note that here we do not consider the contri- 
butions from the operators 05,7,8! they are not only color suppressed but also of order 1/mf. 



C. Soft non-factorizable contributions arising from the intrinsic gluon inside the 
B meson 

The marginal result of Fig. Ha. 2) with the b splitting, namely, b — > bg which is governed by the 
GLAP evolution equation if B is boosted to the infinite momentum frame, implies that the extrinsic 
gluon effect is negligible. In this subsection, we shall study the intrinsic gluon contributions to the 
decay amplitudes. The leading diagram is displayed in Fig. HJa). Here the intrinsic gluon plays 
the spectator role so that this contribution cannot be perturbatively calculated. Before proceeding, 
four remarks are in order, (i) The effects for intrinsic gluon being collinear with the spectator 
quark have already been considered in the transition form factor, (ii) The possible effect that the 
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p 



(a) (b) 

FIG. 4: (a) The contribution to the B — * foK decays arising from the intrinsic gluon within the B 
meson and (b) the diagrammatic illustration to the correlation function given in Eq. ()2.32j) . 

gluon content of /o is produced from the spectator quark is included in the present study since the 
spectator quark is a soft object and therefore the gluon could be one of the constituents of the B 
meson, (hi) As illustrated in the previous section, the extrinsic gluon effect is quite small and hence 
there is no double counting problem, (iv) The amplitude of finding an intrinsic gluon within the 
B meson due to the quantum fluctuations are suppressed by A/m&. It should be stressed that the 
idea of twist expansion for wave functions is not suitable in the present case [see Fig. Eta)]. Since 
there is no hard part ready for a perturbative calculation, the argument for the suppression by 
higher twist distribution amplitudes is no longer justified. Although Fig. Et a ) cannot be calculated 
directly, we can turn to Fig. Et D ) where the B meson state is replaced by a current operator with 
large — p 2 flowing through. In this way, Fig. Etb) can be perturbatively calculated and is related to 
Fig. Et a ) via the reduction formula and quark-hadron duality. 

Let us embark on the calculation. Consider the four-quark operator 5^7^(1— 75)^/3 9/37^(1^75)90: 
in Fig. EI where a and (5 are color indices. The result of Fig. Et a ) can be represented as 

(70^1^(1-75)^ ^7 M (lT75)fe|^-) Fig .H a) = (hK-\20 8 T iq) \B g ) Fig m , (2.30) 

where 

0% q) = 57^(1 - l5 )T a b ff(l T 75)^9, (2.31) 

and we have added a subscript u g" to B~ to emphasize the intrinsic gluon for this case. This 
matrix element can be calculated by considering the following correlation function: 

Uf(p,q) = 1 J d 4 xe-^(^(^),/o(p/o)|T(20|(0),i6(x) 5s G^(x)7 Q n(x)|0) 

= U lG (p 2 ,p' 2 ) Pfl + --- , (2.32) 

where the ellipses (and the following one) denote terms with the structures {p+pkY and (p+Pf Q )^, 
respectively, which are not relevant to the present study, and the superscript "IG" stands for 
intrinsic gluon. 
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In the deep Euclidean region of p 2 , as depicted in Fig. Efb), the correlation function can be 
perturbatively calculated in QCD and expressed in terms of LCDAs of the kaon and /o, 

1 



■qIG(QCD) 



e aPm I d^xe-^(fo( P f )\q(Oh U aT^)G a ^x)q(0)\0) 



6 



x{K(p K )\s(Q) 7 „(l- 75 ) 



d Ak ikxiift-mb) 



(2tt)< 



k 2 



mt 



T^fqfK 



du(j)K(u) / Da 



mt 



1 



1 - ag)4> q ± 



ml - (p-up K -a gPfo ) 2 
Pfj,-\ • 



(2.33) 



To estimate (foK\20^ 9 ^\B g ), we apply the quark-hadron duality to the correlation function. Then 
the ground state contribution to the correlation function can be written in the dispersive represen- 
tation as 

Imn IG(QCD) (s0) 



so I m nIG(phys)/ Si0 ) 
as k 



ds 2 v - ' , (2.34) 

% - r «u % S p 

where so is the threshold of higher resonances. On the left hand side of above equation, the spectral 
density Imn IG ( ph y s ) (s, 0) is given by hadronic contributions and reads 

Imn IG(phys) (s,0) = f B 5 2 B {foK\20%\B g )5(s - m%) + higher resonance states, (2.35) 

where (B~ (pB)\ib'y a g s G afl u\0) = Jb^Pb^- On the right hand side of Eq. (|2.34|) . the spectral 
density ImIl IG ^ CD )(s, 0) can be easily obtained from Eq. (|2.33|) . 

Equating Eqs. ()2.32j) and (|2.33|) and performing the Borel transformation yield 



B 



B 



2 9 
m B ~ P 



m 



exp 



B 



M 2 



mt 



(p-up K - a g p fo ) 



a,, 



exp 



1 + ua g m 2 



a g M 2 



(2.36) 



It leads to the light-cone sum rule 

'qK\20% 
fKf q m, 



(f K\2O 8 T {q) \B g ) 



3f B S 2 



1 



B J 



duipxiu) / da g @la 



mt 



s - umi> Jo 



1 > 



doL q —e v ' B 
a g 



'-ml)/M 2 



a n 



(2.37) 



In Fig.0 we plot (f K\2O 8 ^\B) function of the Borel mass squared. The relevant normaliza- 

*2 X 2 



tion scale for the present process is of order \Jm? B — m 2 2.1 GeV. Using the parameters f q , 8j q ,8 B 
given in Eqs. (|JT0|) . B , (lD4l) . f K = 160 MeV, f B = 185 MeV, m b = 4.85 GeV for the pole mass 
of the b quark, and s = 34 ± 1 GeV 2 Q , we obtc 



tain 



{foK\20^ s) \B g ) = (-0.20 ±0.06) cos 0, 
(foK\20^ n) \B g ) = (-0.23 ±0.06) sin^, 



(2.38) 
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(GeV^) 



FIG. 5: (f K\2O^ q) \Bg) as a function of the Borel mass squared. The solid curve stands for 
(f K\2O^ s) \B g ), while the dashed curve for (f K\2O^ n) \B g ). 



in the Borel window 9 GeV 2 < M 2 < 12 GeV 2 Uf 



Note that since (f K\2O 8 _ (q) \B g 



(foK\20+ \B g ), we abbreviate as O 8 ^ in the ensuring study. In summary, the total con- 

tributions arising from the intrinsic gluon in the B meson to the decay amplitudes are (in units of 

G F /V2) 



A(B- -» f K-) Figm 



A(B* -> f K% iM 



-VtbV t 



(C4 + C 6 ) (V2(f K\2O 8 ^\B g ) + (f K\2O 8 ^\B g 



+ }-(c 8 + c 10 )(-^=(f K\2O 8 ^\B g ) - (f K\2O 8 ^\B g 



(2.39) 



III. RESULTS AND DISCUSSIONS 

To proceed the numerical calculations, we shall follow [ill . Il^ | for the choices of the relevant 
parameters needed for QCDF calculations except for the form factors and CKM matrix elements. 



For form factors we shall use those derived in the covariant light-front quark model [23] . For CKM 
matrix elements we use the Wolfenstein parameters A = 0.801, A = 0.2265, p = 0.189 and f\ = 0.358 
[3] . For endpoint divergences encountered in hard spectator and annihilation contributions we take 
the default values pa = Ph = [see Eq. (|2,9|) ], For the running current quark masses we employ 

m b {m b ) = 4.2 GeV, m&(2.1 GeV) = 4.95 GeV, m b {l GeV) = 6.89 GeV, 
m c {m b ) = 1.3 GeV, m c (2.1 GeV) = 1.51 GeV, 

m s (2.1 GeV) = 90 MeV, m s (l GeV) = 119 MeV, (3.1) 

and m q (p)/m s (p) = 0.0413 16]. The strong coupling constants are given by 

a s (2.1 GeV) = 0.293, a s (l GeV) = 0.489 . (3.2) 

It is ready to perform numerical calculations. At the scale p = 2.1 GeV, the numerical results 
for the relevant a^(foK) and ag 8 (K/o) are 

a u A = -0.0533 - £0.0183, a% = -0.0610 - £0.0064, 
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c% = -0.0568 - iO.0163, a L Q = -0.0612 - i0.0039, 

a% = (74.4 - i4.5) x 1(T 5 , a c 8 = (73.6 - i2.3) x 1(T 5 , 

a^ = (561 + i90) x 10~ 5 , a£ = (560 + i92) x 10~ 5 , 

oi = 1.337 + *0.0288, < 8 (#/o) = < 8 (/o^)- (3.3) 

Note that the effective Wilson coefficients a\ , 04 and aio receive large contributions from the hard 
spectator interactions H(JqK). Consequently, 04 and a§ are similar. Using the distribution am- 
plitudes of the kaon and /o(980) given in Eqs. (|2.16|) and (|2.14|) . respectively, the annihilation 
contributions shown in Eq. (|2.8() can be simplified to 



4 (9) « vra s 



^(X A + 32-^ 2 )-4^^Xi 
2 x m b 



2^(Xi - 2X A + ivr 2 ) - Zr«B^{X\ - AX A + 4) 
mb 3 A 



4 (9) PS QTTOi s X A 



> m ^(2X A -l)+r^B{ q \6X A -ll) 



(3-4) 

m fe ^ J 

where the endpoint divergence X A is defined in Eq. (j2,9j) and only the first Gegenbauer polynomial 
in the distribution amplitudes is kept for simplicity. 

Shown in Fig. ©are the branching ratios of B~ -> f (980)K- and B° -» / (980)K° versus the 
strange-nonstrange mixing angle 9 of /o(980). It is evident that the coefficient B^ appearing in 
the distribution amplitude of /q [see Eq. (|2.14|) ] is preferred to be negative so that the annihilation 
terms make constructive contributions. This indicates the reliability of the sum-rule calculation 
of the fo distribution amplitudes (see Appendix B). When 9 = 0, fo is a pure ss state and hence 
the penguin diagram Fig. 1(a) does not contribute (i.e. the form factor F vanishes). On the 
other extreme with 9 = 90°, fo is purely a nn state and the penguin diagram Fig. 1(b) vanishes 
(i.e. f s = 0). Since (a^ — r^a^) is positive, it follows from Eq. (|2,4|) that, for a finite mixing 
angle, the interference between a^foK) and a^(Kfo) penguin terms arising from Figs. 1(a) and 
1(b), respectively, is destructive for < 9 < ir/2 and constructive for ir/2 < 9 < ir. As stated 
before, the fo — cr mixing angle is slightly favored to be in the second quadrant, i.e. tt/2 < 9 < it. 
It is evident from Fig. H3that this mixing angle solution is also preferable by the measurement of 
B f K. More precisely, B(B~ foK~) is obtained to be (5.5 - 8.2) x 10~ 6 for 25° < 9 < 40° 
and (7.8 — 10.9) x 10 -6 for 140° < 9 < 165°. However, even the maximal branching ratio 11.1 x 10~ 6 
occurring at 9 « is still too small by around 40% compared to experiment. It should be remarked 
that our results for B(B — > foK) are larger than the previous calculations [12I ll^ | owing to the 
large scalar decay constants f s>n we have derived (see Appendix A). 

The fact that the observed /o(980)-ftT rate is higher than the naive model prediction calls for 
some mechanisms beyond the conventional short-distance approach. Some possibilities are: 

• Final state interactions. The predicted B — ► ttK rates in the short-distance approach are 
in general smaller than the data by around 30% (see e.g. ^ ^ s a ^ so known that 

the QCDF predictions for penguin-dominated modes such as B — > K*tt, K<p, K*<p are 
consistently lower than the data by a factor of 2 to 3 Long-distance rescattering via 

charm intermediate states (or the so-called charming penguins) will not only enhance the 
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_l ^- ..J-' . >-J jXT I I I I i_ 



180° 



(a) 



(b) 



FIG. 6: Branching ratios of B — > fo(980)K versus the mixing angle 9 of strange and nonstrange 
components of /o(980), where the solid curve corresponds to the Gegenbauer moments B n ^ s given 
in Eq. (|B12|) and the dashed curve is for the same Gegenbauer moments but with opposite signs. 
Calculations are based on QCD factorization. 



aforementioned penguin-dominated decays but also drive sizable direct CP violation observed 



recently in B° — > K + ir~ mode 48]. The same rescattering effects may be expected to enhance 
/o(980)if rates. Unfortunately, we are not able to estimate the long-distance rescattering 
contributions to the f$Ks rate from intermediate charm states owing to the absence of 
information on foDD and /qD^D*^ couplings. 

Gluonic coupling of the scalar meson. It is known that a possible explanation of the enormous 
production of B — > r(K and B — > vj 'X s may be ascribed to the process b — > s + g + g and 
the two gluons fragment into r/ The same mechanism may be also responsible for the 
enhancement of / (980)K 0- 



Large weak annihilation contributions. Just like the pQCD approach 51] where the annihi- 
lation topology plays an essential role for producing sizable strong phases and for explaining 
the penguin-dominated B — > VP modes, it has been suggested in 0] (see also 0]) that 
a favorable scenario (denoted as S4) for accommodating the observed penguin-dominated 
B — ► PV decays and the measured sign of direct CP asymmetry in B° — » K~ir + is to have 
a large annihilation contribution by choosing pa = 1, <j>A = —55° for PP, 4>a = —20° for 
PV and (j)A = —70° for VP modes. The sign of 4>a is chosen so that the direct CP violation 
A K - n + agrees with the data. Using the same set of pa and 4>a as the PP modes, we found 
B(B~ — ► foK~) = 17.5 x 1CU 6 in agreement with the data. However, the origin of these 
phases is unknown and their signs are not predicted. Moreover, since both annihilation and 
hard spectator scattering encounter endpoint divergences, there is no reason that soft gluon 
effects will affect only pa but not ph- 

Subleading corrections arising from the three-parton Fock states of the fo and/or from the 
intrinsic gluon inside the B meson. It has been shown that the subleading corrections to 
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25 




FIG. 7: Branching ratios of B — » /o(980)-ftT versus the mixing angle 9 of strange and nonstrange 
components of /o(980), where the solid curves are for full decay amplitudes, dashed curves for short- 
distance QCD factorization amplitudes, and only the effects from the intrinsic gluon within the B 
meson are evaluated in the dash-dotted curve. The bands correspond to the errors in Eq. (|2.38j) . 



QCD factorization due to the extrinsic softer gluon effect could enhance the branching ratio 



of Krf to the level above 50 x 1(T 6 [43]. 



In the present work, we have examined in Sees. II B and II C the subleading effects mentioned 
in the last item. For corrections due to the three-parton Fock states of the fo or kaon, we find 
that they are of order and numerically negligible in B — > fo(980)K owing to the fact that 

(/o (980) | Vn\0) = 0. The detailed calculations can be found in Sec. Ill Bl and are summarized in 
Eq. (|2.29|) . As for the corrections originating from the intrinsic gluon within the B meson, they 
belong to the higher Fock components arising from quantum fluctuations and are suppressed by 
A/mfe. Since this mechanism is nonperturbative, the twist expansion for the fo (or K) distribution 
amplitudes is no longer suitable in this case. Therefore, the contribution shown in Fig. Ufa) is 
expected to be suppressed by \/{mi ) ) n with 1 < n < 2. Based on the reduction formula and quark- 
hadron duality, we are able to estimate the nonperturbative effects in Fig. 0(a). A detailed study 
of intrinsic gluon effects has been shown in Sec. Ill Cl and is summarized in Eq. (|2.39|) . 

Taking into account the aforementioned subleading corrections, we show in Fig.[7|the branching 
ratios of B — > fo(980)K versus the mixing angle 9 of strange and nonstrange components of /o(980). 
At 9 ~ 20°, the intrinsic gluon effects can make 25-40% corrections to the decay amplitude so that 
the resulting branching ratio is of order (12 ~ 20) x 10 -6 , in agreement with experiment. These 
corrections are constructive for 0° ^$ 9 < 80° as well as 160° ^ 9 ^$ 180° and destructive for 
80° ;$ 9 ^$ 150°. In short, when the effects due to the intrinsic spectator gluon within the B meson 
are inlcuded, 0° < 9 < 55° and 155° < 9 < 180° are preferable by the B -> f K data. 
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IV. CONCLUSIONS 



In the present work we have studied the decay B — ► /o(980)-K" within the framework of QCD 
factorization. Our conclusions are as follows: 

1. While it is widely believed that /o(980) is predominately a four-quark state, in practice it is 
difficult to make quantitative predictions onB^ foK based on the four-quark picture for 
/o(980) as it involves not only the fo form factors and decays constants that are beyond the 
conventional quark model but also additional nonfactorizable contributions that are difficult 
to estimate. Hence, we shall assume the two-quark scenario for /o(980). 

2. There are two distinct penguin contributions to B — > foK and their interference depends on 
the unknown mixing angle 8 of strange and nonstrange quark contents of /o(980): destructive 
for < 8 < it/ 2 and constructive for ir/2 < 8 < it. A mixing angle in the second quadrant is 
preferable by the measurement of B — > fo(980)K. 

3. Based on the QCD sum rule method, we have derived the leading-twist light-cone distribution 
amplitudes of /o(980) and the scalar decay constant f q . It is found that f s is much larger 
than the previous estimate owing to its scale dependence and the large radiative corrections 
to the quark loops in the OPE series. The measured B — ► JqK rates clearly favor the sign 
of the fo distribution amplitudes predicted by the sum rule. 

4. Based on the QCD factorization approach, we obtain B(B~ — > foK~) = (5.5 — 8.2) x 10~ 6 
for 25° < 8 < 40° and (7.8 - 10.9) x 1(T 6 for 140° < 8 < 165°. Hence, the short-distance 
contributions are not adequate to explain the observed large rates of foK~ and foK° . 

5. Possible subleading corrections from the three-parton Fock states of the fo and from the 
intrinsic (spectator) gluon inside the B meson are estimated. It is shown that while the 
extrinsic gluon contribution to B — > foK is negligible, the intrinsic gluon within the B 
meson may play an eminent role for the enhancement of fo(980)K. 
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APPENDIX A: DETERMINATION OF THE SCALAR COUPLING OF f 

To determine the scalar decay constant f q of /o(980) defined by (0|gg|/o) = f q with /q 
nn = (uu + dd)j\pl and /g = ss, we consider the two-point correlation function 



n(^) = i / d i xe iqx (0\T(f(x)j^(0)\0) , (Al) 

where j g = qq. The above correlation function can be calculated from the hadron and quark- 
gluon dynamical points of view, respectively. Therefore, the correlation function arising from the 
lowest-lying meson /q can be approximately written as 

( / )2 f 2 1 r s o TttiTT ope 

1 1 ds 1 ^-^, (A2) 



where n OPE is the QCD operator-product-expansion (OPE) result at the quark-gluon level, so is 
the threshold of the higher resonant states and the contributions originating from higher resonances 
are approximated by 

- / ds =-. A3 

71" J s s- q l 

We apply the Borel transformation to both sides of Eq. (|A2|) to improve the convergence of the 
OPE series and suppress the contributions from higher resonances. Consequently, the sum rule 
with OPE series up to dimension 6 and 0(a s ) corrections reads 



m { fiy< 2/M2 ( ^#V /b = 

■"' ' \a s (M) j 8vr 2 



a s (M) (n 1(1) M 2 ' 



+ + 3{mqq - q) _ ^± (^q a ^T a qqa^T a q) + ^(q^qq^T^ , (A4) 

where /(l) = 1 - e~ s °/ M2 (l + s /M 2 ), 1(1) = J^ SQ/M 2(lnt)ln(-lnt)dt, and we have taken into 
account the scale dependence of f q , 

4/6 

a s (M) 



f q (M) = f g (J \ . (AT) 



with b = (lliV c — 2n/)/3. Here, the anomalous dimensions of a s G 2 and m q qq are equal to zero, 
while the anomalous dimensions of the 4-quark operators have been neglected. In the numerical 
analysis, we shall use the following values for vacuum condensates and quark masses at the scale 
of fi = l GeV 0: 

{a s G% v GW) = 0.474 GeV 4 /(4vr) , 

{uu) = {dd) = -(0.24 GeV) 3 , (ss) = 0.8(uu) , (A6) 

(m u + m d ) /2 = 5 MeV , m s = 119 MeV , 

and adopt the vacuum saturation approximation for describing the four-quark condensates, i.e., 

(OlqTiTOqqTiTOqlO) = -_L_Tr(r i r i )Tr(T a T a )(^) 2 . (A7) 
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FIG. 8: and f q as functions of the Borel mass squared M 2 . The solid curve is obtained for 

j s = ss and the dashed curve for j n = nn. 



Taking the logarithm of both sides of Eq. (|A4|) and then applying the differential operator 
M 4 d/dM 2 to them, we obtain the mass sum rule for flf\ In Fig. |H1 we explore two scenarios 
for (i) qq = ss and (ii) qq = nn. Numerically, we get 



mjj ~ (1.02 ± 0.05) GeV for s = 2.6 GeV 2 , 



m 



(n) 
fo 



~ (0.99 ± 0.05) GeV for s = 2.6 GeV 2 



(A8) 



where the value of so is determined when the maximum stability for the mass sum rule is reached. 
Substituting the above results for the threshold sq and the masses obtained in the mass sum rule 
into Eq. (|A4[) . we obtain the sum rule for f q at \x = 1 GeV as a function of the Borel mass squared 
M 2 (see Fig. |KJ) . The results are 



and 



f a (ji = 1 GeV) = 0.33 GeV , f n (u = 1 GeV) ~ 0.35 GeV , 



f a (li = 2.1 GeV) ~ 0.39 GeV , f n (u = 2.1 GeV) ~ 0.41 GeV . 



(A9) 



(A10) 



Evidently, they are much larger than the typical decay constant of pseudoscalar mesons^ 
Two remarks are in order. First, the scale dependence of f q was not considered in 



37 



381. If 



the scale dependence is neglected, we will obtain sq ~ 1.6 GeV 2 , in accordance with the results in 
37, 



38j. Second, the large radiative correction to the quark loop in the OPE series arises mainly 



from one-gluon exchange and is likely the effect of the color Coulomb interaction, as in the case 
for the B meson 54] . The values of f q and so , rather than , are very sensitive to the radiative 
corrections to the quark loop and to the scale dependence of f q . However, these effects were not 
considered in 



37 



. If neglecting ce s corrections to the mass sum rule while keeping sq = 2.6 GeV 2 , 
the resulting f s (fj,=l GeV) will be 0.26 GeV. Moreover, if further setting sq = 1.6 GeV 2 , the sum 
rule yields / s (/x=l GeV)=0.19 GeV, in agreement with |37J, |38(]. The validity of our results is thus 
realized. It is interesting to note that a larger scalar decay constant for Kg (1430) has been obtained 



m 



55 



56j based on the sum-rule method. 
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APPENDIX B: LEADING TWIST LCDA FOR f 

The LCDAs & q (x,fi) corresponding to the ideal states /q = qq are defined by 

(/ 9 (P)I9(^)7^(0)|0) = pj q f 1 dxe™P- z <f> g (x,»), (Bl) 

Jo 

where x (or x = 1 — x) is the /o momentum fraction carried by the quark q (or antiquark q) and 
fi is the normalization scale of the LCDA. §> q (x,(j,) can be expanded in a series of Gegenbauer 
l!0 



polynomials 



<&,(*, M) = 6x(l - x) X! ^(M)Cr (2x-l), (B2) 
2=1,3,5,... 

with multiplicatively renormalizable coefficients (or the so-called Gegenbauer moments): 

Consider the following two-point correlation function 

Ik(q) = i J d^xe^iOmOtix) O\0)\0) = (zq) 1 * 1 ^), (B4) 

where 

(0|O z |/ 9 (p)) = (0\q j{iz D) l q\tt(p)) = {zp)J q f (2x - l) 2 %(x)dx = {zp)fM) , (B5) 





<0|O|#(p)> = (0\qq\fg(p)) = mfU (B6) 

with z 2 = 0. The Gegenbauer moments c/)j(x,p,) can be easily expressed in terms of the above 
defined moments (£ l q ), for which the sum rule reads 

- ^ <IM ' + 11 (<*> + + ' < B7 » 

The conformal invariance in QCD exhibits that partial waves, in the expansion of $ g (x,/x) in 
Eq. (|B2|) . with different conformal spin cannot mix under renormalization to the leading-order 
accuracy. As a consequence, the Gegenbauer moments (pi in (|B2|) renormalize multiplicatively: 

^) = U^)(^4V il)l \ (B8) 



where the one-loop anomalous dimensions are 



58 



a s (/x) 



^^'('-(iTW^^gj)' (B9) 

with CV = (AT? — 1)/(2N C ). We consider the renormalization-improved sum rule of Eq. (|B7|) . where 



the anomalous dimensions of relevant operators can be found in 4J], such that 

{qq)n = (qq)l 



4 



(fi^cr • = (g s qa ■ Gq)^ 

2/-y2\ / 2 r~i2 



__2_ 

Qs(Mo) \ 36 



(giG\ = (giG') lt0 . (BIO) 
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FIG. 9: The leading twist-2 LCDAs for the /o meson vs the momentum fraction carried by the 
quark, where the solid curves are evaluated with considering the Bi,B%, B§ central values given in 
Eq. (|B12j) . while the dash-dotted and dashed lines are obtained by setting B§ = and B% = B$ = 0, 
respectively. The darker lines correspond to the normalization scale \i = 2.1 GeV, while the lighter 
lines to /i = 1 GeV. 



In the numerical analysis, we choose the Borel window 0.9 GeV 2 < M 2 < 1.2 GeV 2 , consistent 
with the previous case, where the contributions originating from higher resonances and the highest 



OPE terms are well under control. We also use the input parameters given in Eq. (|A6|) and 44] 

{g s uaGu) = (g s daGd) = — 0.8(mt), 

(g s saGs) = 0.8(g s uaGu). (Bll) 

It should be noted that in the large I limit for the moment sum rule, the actual expansion 
parameter is M 2 /l. As a result, for I > 7 and fixed M 2 , the OPE series are convergent slowly or 
even divergent, i.e. the resulting sum-rule result becomes less reliable. We thus obtain the first 
three non-zero moments, at the normalization scale \x = 1 GeV, 

= -0.55 ± 0.05 => B[ n) = </>" = 5(£*)/3 = -0.92 ± 0.08, 
(f n ) = -0.43 ± 0.02 Bi n) = $ = ^(f n ) - = -1.00 ± 0.05, 

(£) = -0.33 ± 0.01 B&> = 4>l = |(7(£) ~ 3<£» = -0-40 ± 0.05, (B12) 

where the sub-(super-)script n denotes fin. Rescaling to the normalization point [i — 2.1 GeV, we 
find 

(f i) = -0.44 ± 0.04 B[ n) = 0i = -0.73 ± 0.07, 
(£) = -0.33 ± 0.02 => B^ n) = (p3= -0.74 ± 0.04, 

= -0.24 ± 0.01 #J n) = <$>l = -0.017 ± 0.005. (B13) 

As for qq = ss, the analysis yields (^> 3 ' 5 ) ~ ((ss) / (uu)){& 3 > 5 ) ~ 0.8(^ 3 ' 5 ) and b['£ >s ~ 0.SB^ B . 
The leading-twist LCDAs 3> s (x, /•*) and ^(^j A 1 ) f° r the /o meson, which vanish in the large /U limit, 
are illustrated in Fig. 
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APPENDIX C: DECAY CONSTANT 5 2 q 

Jo 

The parameter 5 2 fq appearing in Eq. (12.231) can be denned through the matrix element of the 
local current operator j% = qj a g s G a /3q, 

(o@(o)|/o 9 b)> = i/^. (ci) 

To calculate the 5f parameter, we consider the following non-diagonal two-point correlation func- 
tion 

flp{q 2 ) = i J d 4 xe^(0|T@(x)^t (0)) | 0) = qp u(q 2 ), (C2) 

where j q = qq as in Eq. ()A1() . In the deep Euclidean region of q 2 , the above correlation function 
can be perturbatively calculated in QCD and the result is given by 

n QCD (? 2 ) = -^(^)(ln(W) + ^) -^(qgsa-Gq). (C3) 



After approximating the higher resonances as given in Eq. I)A2|) and performing the Borel trans- 
formation, the resultant sum rule with OPE series up to the order of dimension 7 and radiative 
corrections of order 0(a s ) reads 



32 . .4 
6 



f « 6 f q m f° [a s (M) [a s (M) / 



where the anomalous dimensions of some operators have been listed in Eq. (|B10j) and that of 
qi a 9sG a pq is -32/(96) 0- 

In the Borel window 0.9 GeV 2 < M 2 < 1.2 GeV 2 where the contributions originating from 
higher resonances and the highest OPE term are well under control, using input parameters in 
Eqs. (HSJ) and (|Bllf) . we obtain 

5}s(fi = 1 GeV) = (0.09 ± 0.01) GeV 2 , 

8} s {n = 1 GeV) = (0.10 ± 0.01) GeV 2 , (C5) 
and, as rescaled to u = 2.1 GeV, 

S 2 f s(u = 2.1 GeV) = (0.08 ± 0.01) GeV 2 , 

6 2 r {fi = 2.1 GeV) = (0.09 ± 0.01) GeV 2 , (C6) 

" 

for j q = ss and fin. The results are depicted in Fig. EH 
APPENDIX D: DETERMINATION OF <5| 

The parameter related the intrinsic gluon content of the B meson is defined by 

(0\u 7 a g s G ai ,b\B) = if B 8 2 B p^ (Dl) 
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FIG. 10: 5 2 f q (u = 1 GeV) as a function of the Borel mass squared M 2 with the same notation as 
Jo 

Fig.0 We have used mjf = 1.02 GeV and J s = 0.33 GeV for the solid curve, while mjj = 0.98 GeV 
and f n = 0.35 GeV for the dashed curve. 

To estimate the parameter 8g, we consider the non-diagonal two-point propagator 

IlJ = ij d 4 xe illx (0\T(b(x)g s G afl (xh a U (x) u(0)tb6(0)|0). (D2) 

Following the same line as in Appendix O we arrive at the 5% sum rule 

2 m b + m u m \~"i . „ , 

6 B- 4m 2 f2 e^^(qg s a-Gq), (D3) 

where the contribution arising from a s {qq) is much smaller than the quark-gluon mixed condensate 
shown in Eq. (|G4j) and hence can be neglected. Using the value for quark-gluon mixing condensate 
in Eq. (|B11|) . fs = 185 MeV, and wij = 4.85 GeV for the pole mass of the b quark, we obtain 

6% = 0.022 ± 0.002 GeV 2 , (D4) 

within the Borel window 9 GeV 2 < M 2 < 12 GeV 2 0. Note that the scale dependence of 5 2 B is 
weak. 
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